Mach-Zehnder interferometry with interacting trapped Bose-Einstein condensates 
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We theoretically analyze a Mach-Zehnder interferometer with trapped condensates, and find that 
it is surprisingly stable against the nonlinearity induced by inter-particle interactions. The phase 
sensitivity, which we study for number squeezed input states, can overcome the shot noise limit and 
be increased up to the Heisenberg limit provided that a Bayesian or Maximum-Likelihood phase 
estimation strategy is used. We finally demonstrate robustness of the Mach-Zehnder interferometer 
in presence of interactions against condensate oscillations and a realistic atom counting error. 
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Atom interferometry [1] with trapped Bose-Einstein 
condensates (BECs) is a very promising tool for the most 
precise measurements. The non-linearity of BECs makes 
it possible to create highly squeezed states, which should 
allow to surpass the classical shot noise limit for the phase 
sensitivity AO = 1/y/N by a factor of up to the 
Heisenberg limit (HL) AO = 1/N [2, 3], where N is the 
number of atoms in the condensates. 

Both Atom chips [4] and dipole traps [5] allow for 
versatile control of trapped BECs, and coherent split- 
ting and interference [6, 7] have been demonstrated. 
The preparation of moderately number squeezed states 
through splitting of a condensate by transforming a har- 
monic potential well into a double-well [8] has been re- 
cently achieved experimentally [9, 10], and it has been 
suggested to use optimal control strategies to create 
highly squeezed states at short time scales [11] exploiting 
the atom- atom interactions. 

However, according to the current literature it is gen- 
erally believed that interactions are detrimental for inter- 
ferometry as they induce phase diffusion [12], thereby de- 
creasing the phase coherence [9, 10, 13-16]. The proposed 
standard solution is making the interactions small by em- 
ploying Feshbach resonances [17, 18] or using state selec- 
tive potentials for internal degrees of freedom [15, 16]. 
This is not always possible, and in many cases not de- 
sirable, because Feshbach tuning requires field sensitive 
states which are on the other hand not ideal for precision 
interferometry. Moreover, residual interactions might 
still decrease the sensitivity. 

In this paper, we analyze the Mach-Zehnder (MZ) in- 
terferometer for BECs trapped in a double- well potential 
in presence of atom-atom interactions. We show that 
the sensitivity is not substantially degraded by the in- 
teractions and Heisenberg scaling can be achieved with 
the resources of number squeezed input states and atom- 
number measurements as readout. Our scheme is ro- 
bust against mechanical excitations of the BEC and finite 
atom number detection efficiency. 
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FIG. 1. (Color online) Mach-Zehnder interferometer sequence 
for a finite phase (a) in absence and (b-c) in presence of 
interactions, visualized on the Bloch sphere (T t — T e — 1). 
(a) A number squeezed initial state (small width along z- 
axis, number squeezing factor £n = 0.2) is transformed into 
a phase squeezed one (small width along equator) by a BS 
(rotation around x-axis). Next a phase is accumulated due 
to an external potential (rotation around z-axis). A second 
BS transforms the state such that the phase is mapped onto 
a number difference, (b) Even for very small interactions 
UoN — 0.1 the number squeezing in the final state is lost, 
(c) For larger interactions UoN = 1, the initial state (here 
a Fock state) gets strongly distorted and winds around the 
Bloch sphere. 



The initial state of the interferometer sequence consists 
of two uncoupled, stationary BECs with number fluctu- 
ations AJ Z [19]. We resort here to a generic description 
characterized by two parameters tunnel coupling ft and 
interaction energy Uo [8], and discuss a realistic model at 
the end of the paper. We first introduce the ideal (i.e., 
non-interacting) MZ interferometer as discussed in [20- 
22]. It consists of two cold atom beam splitters (BS) with 
Hamiltonian H t — —VtJ x ^ and in between a phase accu- 
mulation due to an energy difference AE between the two 
wells (with Hamiltonian H e = —AEJ Z ). We visualize a 
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typical interferometer sequence on the Bloch sphere [13] 
in Fig. 1 (a). A BS corresponds to a 7r/2-rotation around 
the x-axis during a time T t . The first BS transforms 
the number squeezed input state into a phase squeezed 
one. The second BS transforms an accumulated phase 
= AET e (a z-rotation caused by an external potential 
during a time T e ) into a relative atom number difference. 
The whole interferometer transformation can be written 
as I^ol^t) = e~ l9 ^ y |^ijv)- A number squeezed input state 
[with squeezing factor £/v := AJ z /(y/N/2) < 1] reduces 
the measurement uncertainty in the atom number of the 
final state [20, 23]. 

Atom-atom interactions are described by the Hamil- 
tonian Hi = UqJ^ [24], and the whole interferometer 
transformation reads now 
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Even for very small interactions, UqN = 0.1, the state 
gets distorted [Fig. 1 (b)]. For larger interactions [Fig. 1 
(c)] the state covers almost the whole sphere. If we em- 
ploy the usual parameter estimation based on the mean 
value of all the measurement results [2, 20, 23], the phase 
sensitivity is worse than shot noise. 

Contrary to the expectations of this estimation, we 
will now show that interactions do not substantially limit 
interferometry. In a completely general fashion we use 
the Quantum Fisher information Fq(\^in)) = 4(Ai?) 2 , 
where R is the generator of the interferometer transfor- 
mation [22, 25], to compute the Cramer-Rao lower bound 
(CRLB), which determines the best possible phase sensi- 
tivity independent of the choice of the measured observ- 
able [26]. For the interferometer transformation Eq. (1) 
we find 
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i.e., it is given by the number fluctuations after the first 
nonlinear BS [27]. m denotes the number of independent 
measurements. 

We start by analyzing a Fock input state |^in) 

{a^ R ) N/2 (a{) N/2 \0). From the scaling of H = H t + Hi 
with N we find AJ z (t = T t ) « aN (with constant a). 
Thus, we expect Heisenberg scaling oc 1/N whenever 
UoN is constant for increasing N [28]. 

Now we have to clarify whether one can indeed achieve 
a sensitivity close to the Heisenberg limit if one is re- 
stricted to a number measurement as in experiments. 
The Classical Fisher information (CFI) [25] 
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allows to estimate a lower bound of AO = 1/ y/mF(0) for 
this specific type of measurement (we consider ^ « 1). 

Hereby, P(n\6) = K^I^out)! 2 * s ^ ne conditional prob- 
ability that an atom number difference n is measured 



for phase 6. Below we choose a constant UqN = 1 and 
vary the BS and accumulation times T t and T e (with 
QT t — 7r/2 fixed). The influence of larger interactions 
can then be extracted through simple rescaling. 

Heisenberg scaling AO = fj/N persists in presence of 
interactions also for a number measurement, see Fig. 2 
(a). The sensitivity is degraded only by an almost N- 
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FIG. 2. (Color online) (a) Scaling of ^mAO with N for 
T t — 2 (circles), T t — 20 (stars), and finite detection error ±2 
(crosses, for T t = 19) and ±5 (diamonds, T t = 21) (UqN = 1 
and T e fixed), compared to shot noise (dashed line) and the 
Uo = HL yfm~A6 = lAy/m/N (dashed-dotted line), (b) 
The prefactor /3 (obtained from fitting) is shown for opti- 
mized T e < 40 (dark solid line), compared to v^^^crlb • N 
(dash-dotted line) and ^/mAO • N for non-optimized values of 
T e < 40 (bright solid line). Also results for a finite detection 
error ±2 are shown (dashed line). 

independent prefactor /?, which varies with T t as is shown 
in Fig. 2 (b) (dark solid line). Fast BSs [T t - 1/(U N)] 
give rise to a prefactor of ~ 1, but also for slower BSs we 
can exploit quantum correlations for MZ interferometry, 
which is relevant for relatively large interactions UqN = 
10 [13]. 

The number readout works because interactions trans- 
form the conditional probability distributions P{n\6), 
which have for the ideal MZ and = a single peak with 
width 1, into a complicated pattern with substructures 
of the same width [see Figs. 3 (c,f,i,l)]. These serve as 
the measurement stick and determine the smallest phase 
which can be resolved [22]. The patterns vary with T e , 
such that some of them show up more distinct l/7V-sized 
peaks [blue line in Fig. 3 (c)], maximizing the CFI of 
Eq. (3) better than others (bright red line). 

Thus, the number measurement is not the 'optimal' 
measurement for all values of T e [26]. We compare the 
maximal and minimal prefactor which can be obtained 
by varying T e [bright and dark solid lines in Fig. 2 (b), 
respectively] . The latter lies close to the CRLB (dashed- 
dotted line). Most importantly, there is no upper limit 
to T e which allows, in principle, to accumulate signals for 
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FIG. 3. (Color online) Probabilities P(n\0) = |HV^ T )| 2 
for (a,d,g,j) binomial, (b,e,h,k) moderately number squeezed 
(£ N = 0.2) and (c,f,i,l) Fock states for N = 100, U N = 1 and 
Tt = 1. The ideal case (no interactions) is shown in (d-f), the 
interacting case in (g-i) for T e = 1, and (j-1) for T e = 10. The 
states are also shown on the Bloch sphere, (a-c) P(n\0 = 0). 



a very long time. 

In many experimental situations only input states with 
finite number squeezing < 1 are available. For the 
ideal MZ the sensitivity increases monotonously with 
number squeezing [black line in Fig. 4 (a)], up to the 
HL v^A<9 = lAy/m/N. 
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FIG. 4. (Color online) (a) Phase sensitivity y/mAO for dif- 
ferent interaction strengths versus the initial number fluctu- 
ations £ N (t = 0) for U N = 1 (N = 100). The symbols 
show results for a simulated Bayesian phase estimation (cir- 
cles: T e = 1, diamonds: T e = 10). (b) Number fluctuations 
after the first BS, £tv(£ = Tt), which determine the CRLB 
Eq. (2). 

We start with analyzing the case of long BS times T t 
(bright red lines). We find a transition of the phase sen- 
sitivity as a function of Starting from ^ = 1 (bino- 
mial state), the sensitivity first decreases up to a point, 
say around £at = 0.2 — 0.3. Then it becomes better again 



and finally approaches the HL for very small values of £at- 
Also the CRLB, Eq. 2, which is a strict lower bound to 
A#, shows a transition. The reason is absence of number 
fluctuations after the first BS whenever the input state 
is only moderately number squeezed [bright red line in 
Fig. 4 (b)]. 

For short BS times T t (blue lines) we find a transition 
only for short phase accumulation time T e (blue dashed 
line). In contrast, a longer T e gives a monotonous be- 
haviour (blue dotted line) similar as the CRLB [blue line 
in Fig. 4 (b)]. 

We can get insight into this behaviour from the con- 
ditional probabilities P(n\0) for different input states in 
Fig. 3. For binomial and moderately number squeezed 
states (d,e), they are close to Gaussian shape [black lines 
in (a-b)]. Interactions wash out the structure of the 
squeezed state (h) and increase the variance in the final 
atom number distribution. Thereby the coherent spin 
squeezing of the initial state is decreased to values even 
worse as compared to the more robust binomial initial 
state [14]. For longer T e , interactions induce substruc- 
tures (k). In contrast, for a Fock input state a complex 
pattern emerges even for very small T e (i), whereas a bi- 
nomial input state does not build up any substructure at 
all (g,j). Visualized on the Bloch sphere, P(n\ff) shows an 
interference pattern whenever a state winds around for 
a long enough time such that it becomes a superposition 
of different phase components [Figs. 3 (k,i,l)]. 

In real experiments, AO as calculated from the CFI 
can be obtained by using a Bayesian (or alternatively 
Maximum-likelihood) phase estimation protocol [29]. 
Thereby a series of m measurements is performed, and 
the atom number difference of each measurement is used 
for the phase estimation. We find that such a proto- 
col gives sensitivity in accordance with the more general 
lower bounds as reported by the symbols in Fig. 4 (a) 
(for m = 20). 

The MZ interferometer is robust against shot to shot 
fluctuations in the atom number or nonlinearity [30]. 
A finite atom counting error has the effect of broaden- 
ing the substructures in the probability distributions as 
P(n\9) oc Y,k perr ° r ( n \ k ) p ( k \°)i where P error (n\k) is the 
error probability for measuring k atoms instead of n. In 
Fig. 2 we show that a binomial error probability with 
width <j = 2 gives rise to just another prefactor, because 
a constant detection error is less important for larger N. 
Even for a detection error a = 5 [9, 31], sub-shot noise 
can be found for N > 2000. 

Implementing the interferometer with trapped conden- 
sates, one achieves the BS by lowering the barrier be- 
tween two split condensates, thereby introducing tunnel- 
ing. The full two-mode physics including the spatial dy- 
namics can be accounted for by the Multi-configurational 
time- dependent Hartree for bosons (MCTDHB) method 
[32] , which represents a framework using time-dependent 
mode functions. For a typical trapping geometry on atom 
chips with ujj_ = 2tt • 2 kHz transverse frequency, we find 
that for tunnel pulses on the order of several millisec- 
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FIG. 5. (Color online) AO for realistic control sequences cal- 
culated with MCTDHB for U N = 0.1 and a binomial state 
with Tt = 8, and UoN = 1 and a highly number squeezed 
state with £zv(£ = 0) = 0.05 and T t = 4. The insets show the 
density (upper panel), the optimal control and tunnel cou- 
pling (lower panel) for the squeezed state and N = 100. 



onds, rapid oscillations are induced in the condensates, 
which lead to unwanted excitations [33]. In our earlier 
work [11, 13] we have developed and demonstrated op- 
timal control [34, 35] within MCTDHB. This allows us 



to design controls for fast BS operations without excit- 
ing the condensates, which is achieved by trapping the 
condensates in stationary states after each of the two 
BSs, while at the same time achieving appropriate tun- 
nel pulses. An approximately 7r/4-tunnel pulse [36] is 
achieved for T t = 4 and highly number squeezed input 
states. In Fig. 5 we show results for binomial and num- 
ber squeezed input states with a phase sensitivity close 
to the HL. 

To summarize, we analyzed the phase sensitivity of a 
trapped BEC Mach-Zehnder interferometer in presence 
of interactions. Heisenberg scaling can be achieved for 
an atom number measurement, and there is no upper 
limit to the phase accumulation time. For finitely num- 
ber squeezed input states the phase sensitivity is char- 
acterized by a transition. We demonstrated robustness 
against condensate oscillations and finite detection error, 
and thus our results can be compared to current experi- 
ments. 
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